Frequency Dependence and Dissipation in the Dynamics of Solid Helium 
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Torsional oscillator experiments on solid 4 He show frequency changes which suggest mass decou- 
pling, but the onset is broad and is accompanied by a dissipation peak. We have measured the elastic 
shear modulus over a broad frequency range, from 0.5 Hz to 8 kHz, and observe similar behavior 
- stiffening and a dissipation peak. These features are associated with a dynamical crossover in a 
thermally activated relaxation process in a disordered system rather than a true phase transition. 
If there is a transition in the DC response, e.g. to a supersolid state, it must occur below 55 mK. 
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As the most quantum of solids, helium is a likely can- 
didate for supersolidity but it was not until 2004 that 
experimental evidence for supersolidity appeared in tor- 
sional oscillator (TO) measurements [2[ . The key ob- 
servation was an increase in the TO frequency at tem- 
peratures below 200 mK, suggesting that some of the 
solid helium decoupled from the oscillator, the "non- 
classical rotational inertia" (NCRI) which characterizes 
a supersolid. However, the onset of NCRI in TO exper- 
iments is broad and is always accompanied by a dissi- 
pation peak, features typical of a dynamical crossover in 
a relaxation process. Does the gradual onset represent 
the DC response associated with a phase transition to a 
new state, broadened perhaps by disorder, or is it a fre- 
quency dependent effect, for example due to flow in a vor- 
tex liquid^] or glass- like relaxation Measurements 
over a broad frequency range can distinguish between 
these possibilities and help identify possible relaxation 
mechanisms, but are almost impossible with TOs. 

We have studied the mechanical susceptibility of solid 
helium directly, by measuring the real and imaginary 
parts of its elastic response (shear modulus and dissi- 
pation), analogous to the TO frequency and dissipation. 
This has advantages over a TO: the measurement is di- 
rect, the modulus changes are large and, most impor- 
tantly, the technique is non-resonant so we could mea- 
sure \x over four orders of magnitude in frequency. We 
see a dissipation peak associated with the previously 
reported 0,@| shear stiffening. The peak and the onset of 
stiffening shift to lower temperatures as the frequency de- 
creases, with an activation energy around 0.7 K. A broad 
range of activation energies is needed to fit the data - the 
relaxation must occur in a highly disordered system. 

The frequency of a TO is usually considered to be 
a direct probe of moment of inertia (and thus of any 
mass decoupling) but things are more complicated when 
the helium is solid. Helium affects a TO through the 
"back action" force it exerts on its walls which could 
be purely elastic or might require a more complicated 
modeljJIa, HJ- The period and damping of the TO (real 
and imaginary parts of the angular susceptibility) are 
determined by the magnitude and phase of this force. A 



decrease in the helium's moment of inertia reduces the 
elastic stress at the walls and increases the frequency. 
An increase in its shear modulus would also raise the TO 
frequency by making the composite TO/helium system 
stiffer, thus mimicking mass decoupling, but this effect is 
too small to explain the observed frequency changes fioj]. 
Elastic forces are in phase with the TO drive and af- 
fect only the real part of the response, but any out of 
phase stresses will produce dissipation. Dissipation, e.g. 
the peak which accompanies 

NCRlQSGilEl, will af- 
fect the TO frequency even without mass decoupling 
However, the frequency change predicted for a simple De- 
bye relaxation peak is nearly always much smaller than 
the observed increase, suggesting that most of the ap- 
parent decoupling may still be due to supersolidity |6j. 
More complicated response functions give better fits to 
the TO data[i], but involve additional parameters. Vary- 
ing the measurement frequency would provide a much 
more stringent test of any model of the helium's response 
but this is difficult with a TO. It has been possible to op- 
erate a TO off resonance, but the frequency range was 
very limited (575 Hz +/- 10%)[13j. The clearest evidence 
of frequency dependence comes from a TO designed with 
two modes [12| - at 496 Hz, the NCRI began at lower 
temperatures than at 1173 Hz. 

Elastic measurements 0jB li, HI on solid helium show 
effects which are clearly related to the TO behavior. The 
shear modulus [i increases by as much as 20%, with the 
same dependence on temperature, amplitude and 3 He 
impurity concentration as the NCRI [7]. This confirms 
the unusual nature of solid 4 He, although these measure- 
ments do not directly address the question of whether 
the low temperature state is supersolid. In this paper we 
extend our shear modulus measurements 0, Hj to much 
lower frequencies, below 1 Hz, and we simultaneously 
measure the dissipation. In some cases, we also measured 
the frequency f r and width of an acoustic resonance, ex- 
tending our frequency range to 8 kHz. 

Measurements were made as in ref. 8. Crystals were 
grown from 4 He (0.3 ppm 3 He) over about 1 hour (using 
the blocked capillary technique) , filling a gap of thickness 
D (0.2 to 0.5 mm) between two transducers with area A 
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(1 cm 2 ). An AC voltage V with angular frequency to — 
27r/, was applied to one transducer (with piezoelectric 
coefficient generating a shear displacement Sx — 

di§V at its surface, a strain e = Sx/D in the helium and 
a stress a on the detecting transducer. This produced 
a charge q which was measured as a current I = u>q, 
giving the shear modulus as fi = ~ = ( A ){jy)- For 
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purely elastic deformations, stress is in phase with the 
applied strain but in general the stress lags the strain 
by a phase angle (j) which is related to the dissipation by 
1/Q = tan <fr <fr for small dissipation). The amplitude 
and phase of the current thus give the real and imaginary 
parts of the helium's shear modulus, although electronic 
phase shifts prevent us from measuring absolute values 
of 1/Q. All data in this paper were taken in the linear 
response regime, at strains e < 10 -8 where there is no 
amplitude dependence or hysteresis [?! [l6j . We measured 
the modulus at frequencies between 0.5 and 2500 Hz. 

Figure la compares the shear modulus /i and acoustic 
resonance frequency f r (for a 33.3 bar sample} with the 
torsional oscillator NCRI (for a 65 bar sample [1]). Figure 
lb shows the corresponding dissipation, determined from 
the phase angle between stress and strain, the width of 
the resonance peak[8[ and the TO amplitude. The shear 
modulus (measured at 200 Hz) and the NCRI (at 910 Hz) 
both increase below 200 mK, with very similar tempera- 
ture dependence, and both have dissipation peaks below 
100 mK. The acoustic resonance shows essentially the 
same behavior. The increase in f r is half as large as the 
increase in \i (3.7% vs. 7.4%), as expected since f r scales 
with the transverse sound speed, i.e. as ^fjl. The acous- 
tic resonance dissipation peak is roughly the same size 
as for the modulus but is broader and shifted to higher 
temperature. Figure 2 shows the frequency dependence 
of fi and 1/Q in a 38 bar crystal. Similar behavior was 
observed in other samples [7J , including the crystal of Fig. 
1. The total modulus change is independent of frequency. 
It is larger, ^ ~15%, than that for the crystal of Fig. 
1 but the temperature dependence is similar. A dissi- 
pation peak is centered near the temperature where the 
modulus is changing most rapidly. The peak and the 
modulus change shift to lower temperatures as the fre- 
quency decreases - the behavior expected for a thermally 
activated relaxation process. For a simple Debye process 
with relaxation time r, the modulus /i and dissipation 
1/Q are related to the real and imaginary parts of the 
shear modulus [17j 
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where /i is the "unrelaxed modulus" (lot 3>1) and 
fi — 8/j, is the "relaxed modulus" (wt <1). The strength 
of the relaxation determines the total change — . For a 

E 

thermally activated process, t(E) = T ef where E is the 
activation energy. The crossover from unrelaxed to re- 
laxed modulus occurs at the temperature where lot = 1; 
at this point the dissipation is maximum and 50% of the 
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FIG. 1: Mechanical response in the 33.3 bar (20.4 
cm 3 /mole) 4 He crystal from ref. 7. a) Shear modulus (at 
200 Hz) and acoustic resonance frequency (left axis) and TO 
NCRI (right axis), b) Dissipation corresponding to Fig. la. 
Solid (red) lines are fits described later in the text. The 200 
Hz dissipation and fit are vertically offset for clarity. 



modulus change has occurred. These points (indicated 
by circles in Fig. 2) can be used to determine the tem- 
perature dependence of r. Figure 3 is an Arrhenius plot 
(1 /T vs log lo) of the crossover temperatures marked on 
Fig. 2. We also show the corresponding points (plus a 
point from the acoustic resonance dissipation peak at 8 
kHz) for the 33.3 bar crystal of Fig. 1. Although there 
is some scatter, the relaxation processes are clearly ther- 
mally activated. The data for the 38 and 33.3 bar samples 
have similar slopes, corresponding to activation energies 
E ~ 0.77 K and 0.73 K, respectively, but are offset ver- 
tically, indicating that they have different attempt times 
t (consistent with sample to sample variations 0) We 
can check whether a simple Debye relaxation describes 
the modulus and dissipation in solid 4 He by comparing 
eqns. (1) and (2) to the 200 Hz data of Fig. 1. The 
dashed (blue) line in Fig. 4a is the modulus at 200 Hz, 
calculated from eqn. (1) using the value of the activation 
energy from Fig. 3 (E=0.73 K) and the modulus change 
= 0.074. The value of t (25 ns) is chosen to make 
the modulus change midpoint agree with the data. It is 
clear that the actual crossover is much broader than for 
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FIG. 2: Frequency dependence of a) shear modulus and b) 
dissipation for an hep 4 He crystal at 38 bar (20.1 cm 3 /mole). 
Curves have been vertically shifted for clarity (the modulus 
scale is for the 2000 Hz data). Dots are transition midpoints 
plotted in Fig. 3. 
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FIG. 3: Frequency dependence and thermal activation for the 
crystals of Figs. 1 and 2. Open symbols are the temperatures 
at which 50% of the modulus change has occurred, closed 
symbols are for the dissipation peaks (as in Fig. 2). Lines 
correspond to activation energies of 0.73 and 0.77 K. 



a simple Debye relaxation. We can get a better fit to the 
modulus data with a smaller activation energy (requiring 
a larger attempt time). The dotted (black) curve in Fig. 
4a is for E=0.175 K and t q =65 [as. Despite the good fit, 
it cannot be correct. First, this value of E is inconsistent 
with the frequency dependence shown in Fig. 3. Sec- 
ond, the predicted dissipation peak (dotted black curve 
in Fig. 4b) is too large by a factor of 2.5. This is also 
true for the fit with E=0.73 K and r D =25 ns (dashed 
blue curve) , which in addition is much narrower than the 
measured peak. No choice of E and r„ can resolve this 
discrepancy since eqns. (1) and (2) imply that the total 
modulus change and the height of the dissipation peak 
are directly related by {^) pea k = 57^- 

However, if the relaxation process has a distribution 
of activation energies rather than a single value, then the 
crossover will be broader and the height of the dissipation 
peak will be reduced (ref. 17, Chapter 4). We can fit the 
200 Hz modulus data using a distribution of activation 
energies with characteristic energy A and width W 



n(E) = Be' 



(3) 



where B is the normalization factor. This distribution 
has a long tail at large E and approaches zero for small 
E. The modulus and dissipation are found by multiplying 
n(E) by the contributions from eqns. (1) and (2) and 
integrating over all energies E. 

The solid (red) curve in Fig. 4a shows the modulus us- 
ing this distribution with A = 0.73K and VF=0.45. The 
fit is much better than for a single activation energy of 
0.73 K (dashed blue curve). The predicted dissipation 
peak (solid red curve in Fig. 4b) has approximately the 
right height and width, although there is additional dis- 
sipation above the peak (nearly constant above 0.25 K). 
This distribution also gives the correct frequency depen- 
dence. The solid (red) curves in Fig. 1 show the predicted 
frequency f r and dissipation for the acoustic resonance 
at 8 kHz (and for /x and 1/Q at 200 Hz) using these pa- 
rameters (A = 0.73K, t =9 ns and VF=0.45). The fit to 
f r is good, showing that frequency dependence retains its 
Arrhenius form with the broadened distribution of acti- 
vation energies. The predicted dissipation peak for the 
acoustic resonance has approximately the correct height 
and width but again there is additional dissipation at 
temperatures above the peak. 

In both TO and shear modulus measurements the 
change in the real part of the response is accompanied 
by damping and the onset temperature decreases as the 
frequency is reduced [Hj], but the dissipation peak is sig- 
nificantly smaller than a simple Debye relaxation would 
predict [fj. This behavior is consistent with a thermally 
activated process if there is a broad range of relaxation 
timesQ. In contrast to TO measurements, the wide fre- 
quency range of our measurements allows us to reliably 
determine the characteristic activation energy for the first 
time and, using this value, A ~ 0.7K, to fit the tern- 
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FIG. 4: Fits of thermally activated relaxation process to a) 
the shear modulus and b) the dissipation at 200 Hz for the 
33.1 bar crystal of Fig. 1. Fits are explained in the text. 



perature and frequency dependence of both the modulus 
and dissipation using a distribution of activation ener- 
gies. The onset of stiffening and the broad dissipation 
peak are not signatures of a phase transition near 150 
mK, but rather of a frequency dependent crossover in 
an activated relaxation process. The temperature of this 
crossover continues to decrease with frequency down to 
at least 55 mK (at 0.5 Hz) - if there is a true phase tran- 
sition reflected in, e.g., the DC elastic response, it must 
be at a lower temperature. 

The effects of annealing on shear modulus and TO 
behavior 0, [ll[ show that defects are involved. It is also 
clear that 3 He impurities are important since the onset 
of stiffening 0] (and of TO decoupl ing|l8|) shifts to lower 
temperature with reduced 3 He concentration (in recent 
experiments with extremely low 3 He concentrations [l9j 
the onset of stiffening was at even lower temperature, 
close to 40 mK) . The large shear modulus changes we ob- 
serve are very difficult to understand except from disloca- 
tion motioned, H3 and the role of 3 He atoms is to bind 
to and immobilize dislocations at low temperatures. As 
the temperature is raised, they unbind thermally, allow- 
ing the dislocations to move and relax the applied stress, 
thus reducing the shear modulus. In this picture, the ob- 
served activation energy is the binding energy for a 3 He 
atom on a dislocation; the value A ~ 0.7 K is consistent 



with estimates from previous experiments 20|, 21 1 . The 
value of t —9 ns, is comparable to the inverse of a ther- 
mal frequency, ksT/h (1 ns at 50 mK). The broad range 
of energies needed to fit our data indicates that the relax- 
ation is occurring in highly disordered network of pinned 
dislocations. This system might exhibit glassy dynam- 
ics, as has been suggested based on the TO response^. 
The "blocked annulus" experiment [l|, which implies that 
the decoupling seen in TO experiments is associated with 
long range order, cannot be explained by local dynamics 
of dislocation motion and remains the strongest evidence 
for supersolidity. Our measurements cannot rule out an 
intrinsic phase transition, e.g. one associated with rough- 
ening due to kinks on dislocations !22|, but it would have 
to occur below 55 mK. 

This work was supported by NSERC Canada. 
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